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ABSTRACT 

Using 3-dimensional hydrodynamical simulations of isolated dwarf spheroidal 
galaxies (dSphs), we undertake an analysis of the chemical properties of their in- 
ner regions, identifying the respective roles played by Type la (SNe la) and Type II 
(SNe II) supernovae. The effect of inhomogeneous pollution from SNe la is shown to 
be prominent within two core radii, with the stars forming therein amounting to ^20% 
of the total. These stars are relatively iron-rich and a-element-depleted compared to 
the stars forming in the rest of the galaxy. At odds with the projected stellar velocity 
dispersion radial profile, the actual 3-dimensional one shows a depression in the central 
region, where the most metal-rich (ie. [Fe/H]-rich) stars are partly segregated. This 
naturally results in two different stellar populations, with an anti-correlation between 
[Fc/H] and velocity dispersion, in the same sense as that observed in the Sculptor 
and Fornax dSphs. Because the most iron- rich stars in our model are also the most a- 
depleted, a natural prediction and test of our model is that the same radial segregation 
effects should exist between [a/Fe] and velocity dispersion. 

Key words: hydrodynamics - galaxies: dwarf - galaxies: evolution - galaxies: abun- 
dances - stars: abundances - Local Group. 



1 INTRODUCTION time in these systems, in order to understand the evolution of 

our own Galaxy. With high-resolution spectrographs on 8m- 
Hierarchical structure formation models predict that mas- , , , ., . , , . , , , ., , , r 

r class telescopes, it is now possible to take detailed spectra or 

sive galaxies formed through the continuous accretion of nu- ■ j- • j i + • u j c i • a u L 

to to individual stars m nearby dwart galaxies, thereby shedding 

merous satellites, and that such a process, at an admittedly , , , . , . , . , . c ,. 

' J light on the latters chemical enrichment and star tormation 

lower rate, should be an ongoing one. Simulations of Milky . . , , , ,. c . , 

' ° ° J histories, and consequently, the detailed accounting ot the 

Way (MW)-like galaxies show stellar halos consisting pre- , , , . c , . 

J v ' b b * number a nd and type ol su pernovae that occurred m the 

dominantly of stellar debris from disrupted satellites that pagt (g _ g _ | M atteucci fc ToslfS . For this reason, a rich 

is kinematically and spatially distinct from the population ,., , - , , , , . . , . , . , , , . 

J v r-^ — II 1 literature exists devoted to studying the chemical evolution 

of surviving satellites (e.g. IBullock et alj 120011 ). This is be- 
cause the survival of a satellite is significantly biased toward 
outlying lower-mass systems with less eccentri c orbits that 
have fallen into the gal axy at later times (e.g iMoore et al.l 
l200rj : ISales et all 120071 ). It is not clear if the star forma- 
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2008) have made chemical models of dwarf spheroidals, em- 



tion histories (SFH) and the gas contents of these systems n • • •, , . , . c ■ 1 c 1 

y ' b I — S| 1 ploying primarily the assumption ot a single- zone framework 

are regulated by internal feedback (e.g. | Dekel & Silk || l986| ) _ the interstellar medium (ISM) and stellar distributions 

and/or (more likely) by the tidal interaction with the MW . . • , ■ 1. 

' v are assumed to remain homogeneous m space throughout the 

and by the ram pressure stripping due to the hot halo of the , , ,. . , . , . , ... 

■' 1 ip ^ calculation. As a consequence, spatial mhomogeneities can- 
Galaxy (e.g. | Mayer et al.| | 2006| ). no ^ ^ e ^ a jj eil m ^ accoU nt by these models, and the values 
If these satellites, or dwarf galaxies, contributed to the of the inferred chemical and hydrodynamica i characteristics 
build-up of the MW, it becomes fundamentally important to can Qnly be interpreted as spatial averages . 
reconstruct how the heavy elements have been produced over 
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Three-dimensional hydrodynamical simulations of the 
feedback between ISM and S Ne within d Sph-s ized ob- 
jects have been carrie d out bv | Mori et all lj2002h and by 
IMarcolini et all l|2006j ). iMori et all |2002l) considered pre- 
galactic objects which later merge into galactic halos; in 
this context, SNe II explode within a few tens of millions 
of years, their energy efficiently coupled to the ISM in such 
a manner as to drive an outflow containing a large fraction of 
the ISM and of the freshly produced metals. Moreover, the 
rapid timescale over which this outflow is effected ensures a 
negligible contribution from SNe IaQ On the contrary, the 
SFHs inferred for surviving dSphs are consistent with being 
less intense and more prolonged (spanning several Gyrs, in 
some cases), suggesting that these objects have been able to 
retain their gas for long periods. 

The simulations of Marcolini et all (|2006l ) accommo- 
dated these prolonged SFHs, allowing for the combined ef- 
fects of both SNe II and SNe la. The reason for such an 
extended SF is due to the effectiveness of radiative losses 
which contrast the rapid expulsion of gas by SNe II, and 
leaves it available for SF. As nowadays dSphs have no gas, 
this scenario requires that the ISM is eventually lost by 
the interaction with the Galaxy, probably through the com- 
bined effects of tidal f orces and ram-pres sure stripping (e.g 
IMarcolini et al. 1 120031 ; iMaver et al1l2006h . Support for this 
picture is giv en by the fact that isolated low m ass dSphs such 
as P hoenix (|Young et al.l I2007T ) and Leo T (|de Jong et al.l 
2008) were able to form stars up to 10 Myr ago. A f urther 
hint in the same direction is given by lHaines et al.l (|2007l ) 
who, analysing a sample of ~30,000 galaxies from the Sloan 
Digital Sky Survey, found that the SFHs of dwarf galaxies 
are strongly dependent on their local environment, the frac- 
tion of passively evolving galaxies dropping from ~70 per 
cent in dense environmen ts, to zero in the rarefi ed field. 

While the model by IMarcolini et al.l (|2006l ) represents 
a noticeable improvement from a hydrodynamical perspec- 
tive, when compared to traditional one-zone semi-numerical 
chemical evolution models, it admittedly considered only 
the elements produced by supernovae. This is not a seri- 
ous drawback, in the sense that the elements synthesized 
and distributed by SNe represent a significant fraction of all 
the elements present today (modulo the contribution from 
low- and intermediate-mass single stars), and the results of 
the model can be readily tested against many observational 
datasets. It has been shown that this mo del is consistent 
with many properties of the Draco dwarf ( Marcolini et alj 
2006) and, with minimal assumptions, the chemical prop- 
erties of the globular cluster a;-Centauri (Mar colini et al.l 
12007ft , which is believed to be the remnant of an ancient 
dSph. 

In the present paper, we analyze the results of the model 
by IMarcolini et al.l (|2006T l focusing on the chemical evolu- 
tion of the stars located in the central region of dSphs. 
It has been shown that in several Local Group dSphs, in- 
cluding Sculptor and Fornax, the chemical properties of the 
stars forming within the galactic c ore can differ substantially 
from those situated a t larger radii l|Tolstov et al.| [2003. 2004; 
Battaglia ct al. 2006). In Section [2 we summarise the char- 



1 Indeed, SN e la ar e explicitly neglected within the simulations 
of lMori et all d2002h . 



acteristics of the lMarcolini et alj (|2006h simulation. We then 
discuss our results in Section [3] compare them to observa- 
tional datasets in Section 2] and draw our conclusions in 
Section [5] 



2 GENERAL PICTURE 



IMarcolini et alj l|2006T ) carried out three-dimensional hydro- 
dynamical simulations of isolated dSphs, emphasising the 
different role played by SNe la and SNe II in the chemical 
and dynamical evolution of the system, in the presence of a 
prolonged intermittent SFH. Although the model has been 
employed to study the general characteristics of dSphs, it 
was built initially to explore the origin and evolution of the 
Draco dSph. In this Section we summarise the main char- 
acteristics of the model and its general behaviour. Details 
of its rationale, as well as the choice o f specific parameters, 
can be found in IMarcolini et all l|2006h . 

2.1 The galaxy model 

At first stars are absent, and the galactic potential is given 
by a static dark halo whose density radial profile (truncated 
at rh,t) is given by: 



ph{r) = ph,o 



1 + 



r 



(1) 



The values of the halo parameters are given in Table [T] and 
lead to a halo mass M h = 6.2 x 10 7 M Q . 

Initially we place the ISM in isothermal equilibrium 
within the potential well, with a temperature of Tism ~ 
T vir = 10 4 K. The initial gas mass is Mis m = 0.18Mh, 
accord ing to the baryonic fraction given by ISpergel et al.l 
(2007), which amounts to 1.1 xlO 7 M . The contribution 
of this baryonic mass to the gravity is neglected. 

As the stars form (see next Section) , they are randomly 
distributed in space following a King profile (truncated at 



P*( r ) = P*,o 



-1 -3/2 



(2) 



After 3 Gyr (the duration of our SFH), a total stellar mass 
A'U — 5.6 x 10 5 Mq is formed. Assuming a stellar mass to 
light ratio in the V band A/*/Lv = 2M0/L0, the mass to 
light ratio inside the stellar tidal radius is Ma/Ly.a = 80, 
in good agreement with observations (e.g. lMatedll998l ). 

The complete list of parameters concerning the galac- 
tic model is g i ven in Table [1] We refer the reader to 
IMarcolini et~afl (|2006h for a detailed description of the 
model. 



2.2 SFH and SN explosions 

The SFH in the model is given by a sequence of 50 istan- 
taneous bursts of identical intensity separated by quiescent 
periods of 60 Myr. The forming stars are randomly located 
within the galaxy following the distribution given by Eq. 
The SNe II associated to the freshly formed stars ex- 
plode at a constant rate for a period of 30 Myr (the lifetime 
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Figure 1. [Fe/H] distribution function of the long-lived stars for the "entire dSph volume" (solid line; r* !C = 130 and r*^ = 650 pc), 
and within two (dashed line) and one core radii (dotted line), at three different times: 1, 2, and 3 Gyrs, respectively. 
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Figure 2. Upper panels: abundance ratio [a/Fe] versus [Fe/H] constructed for a random sample of Ng = 500 stars at three different 
times, for the entire dSph simulated volume. Lower panels: as above, but for stars forming within one core radius. Details on the sampling 
method are given in footnote 2. 



of a 8 Mq star, the least massive SN II progenitor) after 
each stellar burst. 

The SN la rate for a single population decreases in 
time after an initial rise. We adopted t he time-dependent 
rate given by iMatteucci fc Rccchi (200l]), according to the 
Single-Degenerate scenario. 

Each SN explosion is stochastically placed into the 
galaxy according to its radial probability proportional to 
the stellar mass within that radius r: P(r) = M*(r)/M*. 

Finally, we assume that each SN II ejects a mean mass 
MsNii.cj = 10 M©, and each SN la ejects M S Nia, e j = 1-4 M Q . 
Every SN II expels 1.0 M of oxygen and 0.07 M of iron 
while each SN la ejects 0.15 Mq of oxygen and 0.74 Mq 



of iron (e.g. iGibson et aTlll997l . and reference therein). The 
explosion energy of each SN of both types is assumed to be 

£sn = 1q51 er g- See 

iMarcolim" et all (20071 ) for more details 
of the chemical evolution implementation. 



2.3 dSphs evolution 

We found that, although the total energy released by the 
SNe II explosions is larger than the binding energy of the 
ISM, efficient radiative losses enable the galaxy to retain 
most of its gas, which thus remains available for the afore- 
mentioned prolonged SFH. 
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Table 1. Galaxy parameters 





rh,c 


fh,t 


Ph,0 j^* r* lC 


r*,t 


p,,o (M/L V )S 


(10 6 M Q ) 


(PC) 


(pc) 


(10" 24 g cm" 3 ) (10 5 Mq) (pc) 


(PC) 


(10 -24 g cm -3 ) 


62 


300 


1222 


4.3 5.6 130 


650 


1.0 80 



a Value calculated inside r#,t. 



The burst of SNe II associated with each stellar burst 
pushes the bulk of the ISM to the outskirt of the galaxy. 
Once the explosions cease (~30 Myr after each star burst 
episode), the ISM flows back towards the centre of the 
galaxy; when the next burst occurs, the gas is pushed out- 
wards again. This oscillatory behaviour leads to a rather ef- 
ficient and homogeneous pollution of the ISM by the SNe II 
ejecta. We note that as the SFH has been fixed a priori, it 
has no direct relation to the gas reservoir within the galaxy. 
However, stars do form during the quiescent periods be- 
tween bursts, when the gas has "settled" , and an a posteriori 
consistency for the SFH is recovered. As the galaxy in the 
present model does not expel the bulk of its IS M, an exter- 
nal c ause, such as the interaction with the MW (|Maver et al.l 
2006), must be invoked at some point during the evolution to 
deprive the galaxy of its ISM (as the observed Local Group 
dSphs are essentially devoid of gas) . 

Given their lower rate, SNe la do not significantly affect 
the general hydrodynamical behaviour of the ISM, but their 
role is relevant for the chemical evolution of the stars. Be- 
cause of their longer evolutionary timescales, SNe la progen- 
itors created in previous starbursts continue to explode dur- 
ing the subsequent quiescent periods, when the gas is flowing 
back into the central region. During these periods the higher 
ambient gas density (together with the lower SNe la explo- 
sion rate) cause the SNe la remnants to be isolated from one 
another, forming chemically inhomogeneous pockets (here- 
after, we refer to these remnants as "SNe la pockets" ) . These 
pockets are "washed out" by successive phases of expansion 
and collapse of the ISM, due to the effects of SNe II, but new 
pockets may form during the quiescent phases between con- 
secutive starbursts. At odd with the ejecta of SNe II, SN la 
debris is rich in iron and deficient in a-elements. Thus, stars 
forming in the SN la pockets possess lower [a/Fe] ratios and 
higher [Fe/H] ratios than those formed elsewhere. This ef- 
fect is particularly important for the chemical evolution of 
the central galactic region, where the SN la rate is greater, 
and higher ambient gas density is achieved. 

In the next section we highlight the differences between 
stars located in the centre and stars located in the outskirts 
of the galaxy, testing our results against several observa- 
tional constraints. 



radii. Initially, when the ISM is rather metal-poor, stochastic 
effects dominate the chemical evolution owing to the still low 
number of SNs explosions and the inhomogeneous dispersal 
of their ejecta; as a consequence, the IDF is rather broad 
and symmetric, as shown in the first panel of the figure. 
With time, more and more metals pollute the ISM and are 
homogeneously dispersed by the cyclic behaviour of the gas 
flow described in the previous section; the newly forming 
stars are more iron-rich, and the IDF peak moves toward 
the high-metallicity side of the distribution. 

Despite the absence of winds, only ~ 20% of the SNe 
ejecta is present within the star forming region. The remain- 
der of the metals, although still gravitationally bound to the 
galaxy, is driven outwards by the effects of SNe explosions 
(mimicking, to some degree, an outflow) and does not con- 
taminate the subsequently forming stars. 

Even with the homogenising action of SNe II, a degree 
of inhomogeneous pollution persists due to the newly created 
SNe la pockets; the iron-rich stars forming in such pockets 
populate the "tail" at high [Fe/H] in the IDF. As apparent 
from Fig. [1] this tail is primarily found within two core radii 
of the galactic centre, and establishes a chemical difference 
between the central and outlying stellar populations. 

As noted in Section [2j the stars forming in the SNe la 
pockets are characterised, in addition to their higher iron 
abundance, by their low values of [a/Fe] (down to ~ —0.5). 
This is apparent in Fig.[2]which shows the distribution in the 
[a/Fe]- [Fe/H] plane of a sample of 500 stars representative 
of the whole stellar population (upper panels) and of the 
stars within one core radius (lower panels) As expected, 
a higher fraction of cv-depleted stars is present in the cen- 
tral region, where SNe la pockets are preferentially found. 
Indeed, the fractions of stars with [a/Fe]<0.1 present after 
1 Gyr within one core radius, two core radii, and the entire 
galaxy are ~20%, ~15% and ~10%, respectively. These val- 
ues decrease with time, and at t — 3 Gyr they are reduced by 
a factor -~30%. This reduction is due to the decreasing rele- 
vance of the contribution of the single pockets, as the mean 
value of [Fe/H] increases with time while the iron content of 
the pockets is always the same. In fact, from Fig. [5] (as well 
as from Fig. [TJ it is clear that the maximum value of [Fe/H] 
never exceeds ~ —0.7. Actually, this value is attained by 



3 RESULTS 

Although our model was tailored to the Draco dwarf galaxy, 
it provides a valuable baseline with which to illustrate sev- 
eral generic properties of dSphs. 

Figure [1] shows the iron distribution function (IDF) of 
the long-lived stars (with masses M < 0.9 Mq) at three 
different times (1, 2 and 3 Gyrs) for the whole dSph (solid 
line) and within two (dashed line) and one (dotted line) core 



2 One of the outputs of our code at any time is the three- 
dimensional stellar distribution tp(f, o, r) (where, for each star, 
/ = [Fc/H], o = [O/Fc], and r is the position). We can ob- 
tain $(/) = J dr J tp(f, o, r)do (shown in Fig. [TJ and <j>(f, o) = 
J ip(f, o,r)dr. The random sampling shown in Fig. [2] is done as 
follows: we first extract randomly a value of [Fe/H] from the dis- 
tribution if), then a value of [a/Fe] is randomly extracted from 
the <j> distribution calculated for the fixed value of [Fe/H]. Similar 
procedures have been followed to obtain the samplings shown in 
Fig. OH |6]and [7] 
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Figure 3. Age-Z, age-[Fe/H] and age-[ct/Fe] distributions for 500 sampled stars within one core radius of the simulated volume. Note 
the metallicity and [o/Fe] spread amongst coeval stars. 



the stars forming within each pocket, where the amount of 
iron delivered by the SN la is constant (0.74 Mg), and the 
amount of H does not vary appreciably because the stars 
form during the quiescent periods, when the ISM settled 
within the galactic core recovers substantially to the same 
densities. As a consequence, the maximum obtainable value 
of [Fe/H] is rather insensitive to the ISM mean metallicity, 
and depends instead on the details of the dSph model. In 
the particular model we are discussing, the SNe la inhomoge- 
neous pollution is expected to become irrelevant at a nearly 
solar value of [Fe/H]. The effects of the inhomogeneity are 
thus expected to be cr ucial in the metal-poor dwarf galaxies 
l|Marcolini et al,ll2007l ). but not in those more metal-rich. 

The stars in Fig. [2] are arranged in two branches, a 
nearly horizontal plateau and an oblique "ramp" connected 
to the plateau by a "knee". The stars located in the upper 
envelope of the plateau are mostly polluted by the ejecta 
of SNe II, while the stars below the plateau, as well those 
on the ramp, are also inhomogeneously (b ut significantly) 
pollu ted in different measure by SNe la (see lMarcolini et al.l 
l200fj , for more details) . As the mean iron content of the ISM 
increases with time, the knee in the [a/Fe]- [Fe/H] diagram 
moves rightward (toward larger values of [Fe/H]) and, to a 
lesser extent, downward (to lower [a/Fe]). 

A further evidence of the inhomogeneous pollution is 
given in Fig. [3] where we plot the a,ge-log(Z /Zq), age-[Fe/H] 
and age- [a/Fe] relations for a randomly-drawn sample of 500 
stars formed within one core radius. In the first two panels a 
striking spread in the metallicity Z and in the iron content 
of the star is apparent, which tends to reduce with time, as 
expected. The [Fe/H] spread is more accentuated because 
stars formed inside SNe la pockets have a much higher iron 
content, but essentially the same Z, than the stars formed 
elsewhere, since SNe la deliver a significant amount of iron 
but a negligible fraction of metals in comparison to SNe II 
(for example, after 3 Gyr the SNe la ejecta accounts for 
30% of the iron produced by all stars, but only for 2% of the 
metals). 
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Figure 4. Abundance ratio [a/Fe] versus [Fc/H] for 500 sampled 
stars of the reference model at t = 2 Gyr, compared to a datasct 
of 28 dSph stars (Draco: triangles ; Fornax: crosses; Leo I: aster- 
isks; Sextans: diamonds; Sculptor: circles; Ursa Minor: squares) 
collected from the literature (see text for details). The stars are 
sampled as in Fig. [2] 



It is commonplace to assume that stars have chemi- 
cal elements in solar proportion, thus leading to an implicit 
equivalence between the IDF and the metal distribution 
function (MDF). That said, a comparison between the first 
two panels of Fig. [5] indicates that t he two distributions are 
different (see iMarcolini et al.l 120071 . for more details about 
this point and how it influences the color-magnitude dia- 
grams) . 

Finally, the third panel of Fig. illustrates the evolu- 
tion of [a/Fe] and shows that the SNe la pockets start to 
form quite early in the evolution (~300 Myr after the first 
star formation episode, as the SNe la rate starts to become 
appreciable) and the tendency of its spread to decrease with 
time (see also the lower panels of Fig. [2} . 
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Figure 5. The radial profile of the 3D velocity dispersion (solid line with crosses) and projected velocity dispersion using fOOO sampled 
stars of our model (thick solid line) obtained solving the Jeans-equation compared to the velocity dispersion of Draco dSph obtained by 
Lokas ct al. I d2005l , left panel) for three sam ples of 207 (dotted), 203 (dashed) and 189 (dash-dotted) stars, and compared to the velocity 
dispersion obtained bv lMuhoz et all (|2005| . right panel). 



4 COMPARISON WITH DATA 

In Fig. [4] we compare the chemical characteristics of the 
stars of the mo del (at 2 Gyr) with data published by 
several authors (IShetrone et alj |200ll 120031 ; iTolstov et all 
120031 ; iGeisler et all 120051 ) for different dSphs (Draco, For- 
nax, Leo I, Sextans, Sculptor and Ursa Minor) observed us- 
ing VLT/UVES. 

Although the dataset is a collection of values of differ- 
ent dSphs, the general pattern of the observed [a/Fe]-[Fe/H] 
diagram is reproduced quite well by our single model. Due 
to the small number statistics in publi shed results it is dif - 
ficult to draw definitive conclusions. ITolstov et ail (|2003l ) 
suggested that the similar abundance patterns across sev- 
eral dSph galaxies was evidence for very similar chemi- 
cal evolution histories, similar initial conditions and simi- 
lar initial mass functions, despite their differen t SFHs, al- 
thou g h there is also evidence for difference s |Venn et al.l 
120041 : ITolstov et alJlioOrj : iLetarte et al.ll2007n . Larger sam- 
ples coming from FLAMES high resolution data sets should 
make this issue clearer (Hill et al. in prep; Venn et al. 
in prep), greatly increasing the statistics inside each sin- 
gle galaxy, including the Sculptor and Fornax dSphs. In- 
deed, one should compare the model with a dataset for 
each dwarf galaxy as the corresponding "knee" can start 
at different [Fe/H] values depending upon the chemical en- 
richment history of the dSph citeptolstoy2006,latarte2007. 
However, for the aim of this paper, the literature data are 
sufficient to highlight the main characteristics (note that the 
future larger samples will overlay the present UVES sample) . 
Here, we simply wish to highlight the point that a-depleted 
stars are observed in several dSphs, just as predicted by the 
model. 



Actually, a-depleted stars can also be explained 
within the framework of th e ch emical model of 
lLanfranchi fc Matteuccil l|2004 120071 ). in which the 
[a/Fe] decrease at larger values of [Fe/H] is the consequence 
of a galactic wind which reduces the amount of available 
gas; the formation of new stars is thus tempered, along 
with the a-elements delivered by SNe II. SNe la, instead, 
having progenitors with l onger evolutionary timescales 
|Matteucci fc Greggiol Il98fj ), continue to inject iron into 
the ISM, lowering the [a/Fe] ratio of the gas. In this case 
the a-depleted stars should be the youngest ones, with 
the largest values of [Fe/H], and exhibiting a small spread 
in age. Preliminary results of [a/Fe] measurements of 
~ 90 stars in Sculptor and ~ 55 stars in Fornax recently 
published by ITolstov et all (2006) seem to indicate that 
a-depleted stars are found at any [Fe/H]> —1.8 down to 
[Fe/H]~ —0.6 inside each galaxy, in good agreement with 
our model. 

In Fig. [S] the projected radial velocity dispersion pro- 
files of the model ar e plotted and tested agai nst data in- 
ferred for Draco by iMunoz et al l (|2005l ) and lLokas et al.l 
l|2005l ). The velocity dispersion profile is calculated solving 
the Jeans equation and assuming velocity isotropy (/3 = 0). 
Different choices of the (3 parameter can change the results 
dramatically, even if there is no valid motivation, a priori, 
for choosing highly tangential or radial orbits; a s such, /3 = 
appea rs a fair conservative value. For example iMateo et al.l 
l|2007l ) found that for Leo I stars, those situated more to- 
wards the centre appear to have isotropic kinematics, while 
only those located outside a "break radius" of about 500 pc 
exhibit significant velocity anisotropy. These authors also 
suggest that the break radius represents the location of the 
tidal radius of Leo I. 
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Figure 6. Left panels: [Fc/H] projected radial distributio ns of 300 s ampled stars at three different times (1, 2 and 3 G yr), Right panels: 
observe d [Fe/H] radial distrib utions for 87 stars in Draco llFaria et alj|2007h . 308 stars in Sculptor llTolstov et al, 2004) and 562 stars in 
Fornax (Battaglia et al. 200a). The empty squares represent stars with a > 8.5 km s _1 while black dots are stars with a < 8.5 km s _1 . 



In the same figure, we also plot the actual 3D ve- 
locity dispersion profile within the galaxy (solid line with 
crosses). While the projected profile is essentially flat (within 
1 km s _1 ) up to 500 pc, the the actual profile shows a central 
minimum, ~ 3km s~ below the maximum value attained in 
the outer regions. The tendency of the central stars of the 
model to have somewhat lower values of the velocity dis- 
persion is emphasised in the left panels of Fig. [6] where we 
show the projected radial positions and the ratios [Fe/H] of 
a sample of 300 stars: here the empty squares represent stars 
with a > 8.5 km s _1 while black dots are stars with a < 8.5 
km s _1 . These diagrams also show that [Fe/H] rich stars 
are found preferentially in the central region of the system 
(and only with lower velocity dispersion), while the main 
Fe-poor population is more smoothly distributed. This re- 
sult can be easily understood: while the iron-rich stars form 
preferentially in the SNe la pockets which are more con- 
centrated toward the center, the stars with lower values of 
[Fe/H] form throughout the galaxy from an ISM whose pol- 
lution does not depend strongly on the radius because of 
the dynamically-driven homogenisation (cf. Section [2}. As 



a consequence, two different stellar populations are found 
in our model, with iron-poor stars showing higher velocity 
dispersion and iron-rich (centrally concentrated) stars show- 
ing lower velocity dispersion. A similar anti-correlation be- 
tween the iron content and th e velocity dispersion has been 
actu ally observed in Sculp tor (jTolstov et al.l[2004l ) and For- 
nax l|Battaglia et al.ll2006l ) dSphs (and found in preliminary 
models that we adapted to these galaxies). 

From the left panels in Fig. [6] it is apparent that the 
metallicity gradients occurring in our model are rather small, 
at variance with monolithic dissipative models of elliptical 
galaxy formation. In these models strong gradients form as 
a consequence of existing radial flows carryi ng inward th e 
heavy elements produced in the outer regions |Larson|ll974l ) . 
In our model, instead, because of the small size of dSphs, 
the gas flows powered by the SNe II tend to homogenize the 
ISM, thus inhibiting the formation of chemical gradients (cf. 
section [2|; the present mild [Fe/H] gradient is mainly due 
to SNe la, as discussed in [3] Indeed, being that the oxygen 
(and a-elements in general) is produced for the most part 
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by SNe II, the [O/H] gradient is even shallower than that of 
[Fe/H] l|Marcolini et al.ll2006l) . 

From the left panels in Fig. [6] it is also apparent that 
the small metallicity gradient (0.2-0.4 dex; see Marcolini et 
al. 2006) tends to decrease with time as the amount of iron 
in the ISM increases and the stars forming at larger radii 
acquire increasing values of [Fe/H]. This tendency seems 
to be confirmed b y the stella r radial distribution s of Draco 
llFaria et all 120071) . Sculpto r |Tolstov et all 12004) and For- 
nax ( Battaglia et al.l I2006T ) (right panels in Fig. [6} : these 
dSphs seem to possess s hallower gradie n ts, as their SFH (e.g . 
iBuonanno et all 1 19991 ; iDolphidliool iTolstov et ail l200i 
and reference therein) is more prolonged. While this be- 
haviour can be described only q ualitatively by the reference 
model of iMarcolini et all l|2006f) . further models tailored to 
Sculptor and Fornax are planned, in order provide a more 
appropriate and direct comparison to the data. Neverthe- 
less, our model seems to be in good agreement with the 
"central concentr ation" of the "small [Fe/H] tail" found by 
iFaria et all (|2007f ) for Draco. 



We stress that the presence in the galactic center of 
an iron-rich stellar population could be explained in prin- 
ciple also by oth er models. In the framework depicted by 
iMori et all |2002h the SNe feedback following the initial star 
formation blows out the gas and stops star formation tem- 
porarily. After SNe feedback becomes weaker, more metal- 
rich gas can return and form a new generation of stars dis- 
tributed closer to the center. As pointed out in Section [T] 
however, the model of IMori et all <|2002l ) is not suitable in 
the case of the prolonged SFHs often encountered in dSphs. 
A more metal rich ce ntral population of stars is also ob- 
tained in the model of iKawata et alj (|2006l ) as the result of 
a steep metallicity gradient within a single population, in- 
duced by dissipative collapse of the gas component. This 
is the same mechanism that is suggested to explain the 
metallicity gradient for larger spheroidals, i.e. normal ellip- 
tical galaxies. It is interesting to note that in these galaxies 
the observed gradient of [a/Fe] has a mean value close to 
zero, an d that it does not correlate with galactic proper - 



ties (e.g. ISanchez-Blazquez et"aill2007j; iPipino et all 



2007) 



Thus, within the scheme proposed by IKawata et al.l (2006), 
dSphs are not expected to present a well defined [a/Fe] 
gradient. On the contrary, in the model of IMarcolini et ail 
2006) a-depleted stars are centrally segregated (see Fig. [7} 
as they are correlated with the Fe-rich stars which are pref- 
erentially found in the central region. Thus the observations 
of the [a/Fe] radial profile in dSphs can pot entially select be- 
tween the m odel of [ Marcolin i et all l|2006h and the scheme 
pr oposed by | Kawata et all ([2006). As an aside, the model 
of IMarcolini et al* ( 20061 ) also predicts that the most a- 
depleted stars have lower velocity dispersions, as illustrated 
in Fig. 



We point out that, although the proximity of the MW 
may change the stellar orbits (expecially in the outskirts) of 
a local dSph, a segregation of the a-depleted stars should 
still be present because such stars are located within two 
core radii, a region only slightly affected by the tid al inter- 
action (e.g. lGao et~aill2004l ; TPenarrubia et alj|2008l ). 
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Figure 7. [a/Fe] projected radial distributions of 500 sampled 
stars at three different times (1, 2 and 3 Gyr). The empty squares 
represent stars with a > 8.5 km s — 1 while black dots are stars 
with a < 8.5 km s" 1 . 



5 CONCLUSIONS 

We analysed the out come of the 3D hydrod ynamical model 
of the Draco dSph of lMarcolini et al] l|2006h , and found that 
the results can provide a useful, if qualitative, description of 
the larger population of generic Local Group dSphs. We can 
summarise our major findings as follows: 

i) Evidence is mounting about the effects of inhomo- 
geneou s pollution on th e chemical history of the stars in 
dSphs (|Koch et al.ll200St ). T he effect of the i nhomo geneous 
pollution of SNe la found in IMarcolini et al. (2006) is more 
important in the central region (within two core radii) of 
dSphs. Here the stars forming in the SNe la pockets can 
amount to ~ 20% of the total. 

it) The chemical homogenisation of the ISM by 
SNe II, together with the combined inhomogeneous pol- 
lution by SNe la, naturally accounts for a radial seg- 
regation of Fe-rich stars (with depleted [a/Fe] ratios) 
in the central regions of dSphs. Stars with this abun- 
dance patt ern have been observed i n the center of sev- 
eral dSphs ilShetrone et al.ll200ll. [20031; iBonifacio et al.ll20ol ; 
iMonaco et al.l 120051 : ITolstov et alj|2006l ; iKoch et alj|2008l ). 
As discussed, the most Fe-rich stars in our model are also 
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Q-depleted, and thus the same radial segregation should be 
observed to test and prove our model. 

Hi) At odd with the projected stellar velocity dispersion 
radial profile, the actual one shows a depression in the cen- 
tral region, where the iron-rich stars are partly segregated. 
This naturally entails two different stellar populations with 
an anti-correlation between [Fe/H] and velocity dispersion, 
which has been obs erved in the case of th e Sculptor dSph 
(iTolstov et alj|2004 ) and the Fornax dSph l|Battaglia et all 
l200fj ). 

In a forthcoming paper we will compare new simulations 
tailored to specific dSphs (e.g. Fornax and Sculptor) with the 
properties of hundreds of st ars observed with VLT-FLAMES 
that will be available soon l|Venn fc Hillll2005l ; [Tolstoy et all 
120061 . and references therein). The new models will take into 
account the yields from low- and intermediate-mass asymp- 
totic giant branch stars, as well as the interaction with our 
Galaxy, in order to describe how the dSphs lose their gas. 
For example, an interaction with the MW seems to be par- 
ticularly important in the evolution of the Sagittarius dwarf. 
Preliminary results show that the iron-rich, Q-depleted stars 
in this galaxy are consistent with a SF strongly affected by 
this interaction, as well as with the formation in situ of M54 
as the core of this galaxy. 
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